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1.0  INTRODUCTION 

The  work  reported  under  this  contract  was  motivated  by  the  need 
to  investigate  possible  LF  systems  that  are  capable  of  providing  means 
to  communicate  with  submerged  submarines.  Since  ULF  waves  have  a large 
skin  depth  in  sea  water,  it  seems  natural  to  consider  them  for  the 
stated  purpose,  \ 

The  basic  problem  encountered  in  designing  any  feasible  LF  com- 
munications system  has  three  major  elements: 

• Generation  of  ULF  signals 

• Propagation  from  source  to  submerged  receiver 

• Detection  of  the  signal 

basis  for  the  present  work  was  formed  in  1974  when  TRW  personnel 
performed  preliminary  calculations  of  generation  of  ULF  signals  by 
satellite-launched  ion  beams.  These  results  indicated  that  signal 
levels  of  some  10  my  at  sea  level  could  be  produced  by  beams  of  some  tens 
of  kilowatts,  emitted  at  3 to  4 earth  radii  altitude  at  the  geomagnetic 
equator.  A brief  description  of  this  work  appears  in  Proposal  No.  TRW 
26766. 000R1  (September  1975)  and  will  not  be  repeated  here. 

'The  model  used  in  the  earlier  calculations  was  a "broad"  ion  beam. 
The  term  "broad"  means  that  finite  beam  radius  effects  were  ignored. 

Thus,  such  an  idealization,  which  is  mathematically  convenient,  must 
lead  to  larger  signals  and  higher  efficiency  than  the  realistic  case  of 
a beam  of  finite  radial  extent.  Since  the  broad-beam  calculation  gave 
encouraging  results,  it  was  considered  worthwhile  to  examine  the  finite- 
cross-section,  or  "pencil"  beam. 

The  results  of  this  analysis  of  a finite-radius  beam  of  ions  are 
presented  ,m_Sections  2 through  6 of  this  report.  Section  7 addresses 
the  estimaH  of  spatial  or  radial  spreading  of  the  beam  due  to  electro- 
static effects,  while  Section  8 contains  a discussion  of  current  and 
charge  neutralization.  Modulation  of  the  beam  is  addressed  in  Section 
9,  and  use  of  heavier  ions  (other  than  protons)  is  treated  in  Section  10. 


Section  11  contains  a large  number  of  conclusions  drawn  by  consideration 
of  the  results  in  Sections  5 and  6.  Finally,  Section  12  contains  the 
recommendations  based  on  conclusions  presented  in  Section  11. 
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2.0  BASIC  EQUATIONS  FOR  RADIATED  FIELDS 

Maxwell's  equations  for  a electromagnetic  wave  field  propagating  in 
a magnetized  plasma  can  be  combined  to  produce  a single  wave  equation  for 
the  electric  field  E (r.,t) 

v(v-|)-v2!^2 !/»=-%  (1) 

C P t 

where  J is  the  self-consistent  current  density  which  supports  the  wave  in 

^ ~ B 

the  background  plasma,  while  ^ is  the  externally  supplied  beam  current  density. 
We  introduce  phasor  notation,  and  use  the  convention 

(2)-h  (E  - ivE  ) v = ± 1 
£ = x y 

V Eq  v = 0 (2) 

where  v = 1 (-1)  corresponds  to  the  left-hand  (right-hand)  circular  polar- 
ization with  respect  to  the  DCmagnetic  field  B0,  taken  to  be  in  the  z-direction. 

'Xi  +■ 

The  subscript  v on  other  quantities,  such  as  current  density  Jv,  wave  number 
k^,  etc  follows  the  same  convention  as  defined  in  equation  (2)  above. 

We  next  introduce  the  Fourier-Laplace  transform 
r , r -(pt  + Ik  • r)  . 

Ev  ) =Jd3r Jdte  v Ev  (r,t)  (3) 

of  (3)  to  (1)  yields  the  algebraic  equation 


f 


The  plasma  current  is  / = n ^ q v where  q = ± e and  v is  the  self- 
consistent  wave-induced  velocity,  the  summation  being  made  over  all  charged 
species  present.  The  velocity  v is  determined  by  the  cold  fluid  equations. 


i = m k + k x « 


where  £ = (qHo/mc)j£  is  the  psuedovector  gyrofrequency . The  transformed 
solution  to  (5)  is 


(q/m)  E 


v p + ivfi 


We  define  the  quantity  q which  has  dimensions  of  a wave  number,  by 


2 2 2 V' 

c % ■ p + 2-  -f— 

± p + ivn 


2 2 

where  u>p  = 4irnoe  /m  defines  the  plasma  frequency  for  species  of  mass  m. 
Substitution  of  (6)  and  (7)  into  (4)  yields  the  following  expression  for  E^: 


^'7T7  [Mt  • V ■ ^ Jv  ] 


6 


We  can  eliminate  the  quantity  k • E if  we  multiply  (8)  by  k_v  and  sum 
the  result  over  v to  obtain 


* • £ 


k o* 

-p  p 

— — 1 

* + % 


I — F 

-p  p 


(9) 


Elimination  of  k • E between  (8)  and  (9)  yields 

% % 


E 

v 


- 4np/c 
~1~ 


(kvg  + 


JB) 

v ' 


(10) 


The  Maxwell's  equations  relating  B + E,  when  transformed,  can  be 
manipulated  to  yield  the  corresponding  expressions  for  the  wave  magnetic 
components 


B 

v 


vC 

P 


k E 


o v 


(v  = ±1) 


B 

o 


y^o 


k E 

-p  p 


(v  =0  ) 


(ID 
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3.0  MODEL  OF  THE  ION  BEAM 


We  assume  that  the  ion  beam  is  injected  at  t=o  with  a well  defined 
speeds  VM  (parallel  to  B0)  and  V (perpendicular  to  B0).  The  ions  spiral 
around  B0  with  a Larmor  radius  R = V./fi,  while  their  Lamor  phase  advances 
as  <t>  = Az  where  A = V||/n  is  the  Larmor  wave  number.  At  a time  t > o the 
beam  then  exists  in  the  region  o < z < V | ( t . In  a plane  perpendicular  to 
B0  the  beam  cannot  be  assumed  to  be  sharply  defined  because  space-charge 
and  high-frequency  electrostatic  beam-plasma  instabilities  must  occur. 

These  effects  will  act  to  diffuse  the  beam  across  magnetic  field  lines, 
until  the  beam  density  drops  to  a value  comparable  to  the  background  density. 

To  account  for  this  expected  radial  spread  of  the  guiding  centers  of 
ions  comprising  the  beam,  we  use  a Gaussian  function  centered  on  the  mean 
gyroradius  R and  having  a spread  A.  The  mathematical  model  of  this  beam 
is  then  given  by  the  distribution  function 

nR  2 -1 

fB  (£.»  v,  t)  = — 6 (vz-V||)  6 (Vj-VjJ  6 U-Ax)  • (itA  ) • 

exp  I -(x- Rcos<f.)2/A2  - (y-Rsin<j>)2/A  j • 

H (t)  [H  (z)  - H (z-V ! ! t)  ] (12) 

where  H (t)  is  the  Heaviside  unit  step  distribution,  6 (x)  is  the  Dirac 
delta  distribution,  and  nD  is  the  number  of  beam  ions  per  unit  length  (the 

b 2x 

volumetric  number  density  is  ng/A  ). 

The  beam  current  density  is 

JB  = e/d3v  v f (jr,  v,  t) 


which  yields  the  phasor  quantity 


v 


eVv  e-1vAz  a (*>  t) 


(13) 


where 


G (r,  t)  = (ttA  ) exp  |-  (y-Rcosij))  /A  - (y-Rsi  n<j>)  /A  J H(t)  | H(z)-H(z-V)(  t)J 


and  V = V1(  , V = V in  (13) . 

O II  ±1  x 


(14) 


In  order  to  maintain  the  charge-neutrality  of  the  spacecraft,  an  elec- 
tron current  must  be  emitted  such  that  it  just  cancels  the  beam  ion  current. 

Therefore,  in  what  follows,  we  assume  that  the  total  current  parallel  to  B0 

B r'J 

vanishes , i .e.  J = 0. 

oil 


4.0  HYDROMAGNETIC  APPROXIMATION 

D 

Under  the  assumption  J 0||  = 0 the  parallel  component  of  the  wave 
electric  field  is  just 


4irp  £ 

c2  k2  + q 2 
o 


(15) 


For  low  frequencies  compared  to  the  electron  plasma  frequency, 
that  is  ip  « w , the  quantity  qQ  ^ wpe/c-  Since  low  frequency  waves  have 
very  long  wavelengths,  then  k <<  qQ  for  the  hydromagnetic  modes  (ip-^  n) 
of  interest  to  us.  If  one  compares  the  order  of  magnitudes  of  EQ  and  E+1, 
one  finds  that  Eo  is  smaller  by  the  ratio  k2c2/u>2e<<  1.  Therefore,  we  can 
justify  the  neglect  of  all  terms  involving  Eq  in  our  subsequent  calculations 
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The  polarization  E is  associated  with  longitudinal  electrostatic  waves 

o 

having  virtually  no  magnetic  components  in  cold  plasma. 

In  order  to  further  simplify  the  calculation,  we  snail  concentrate  on 

waves  having  frequencies  much  less  than  the  ion  gyrofrequency . In  practice, 

2 2 

this  is  expressed  by  p « n . Even  |p|  -v  n/2  will  satisfy  the  condition  to 

within  a reasonable  accuracy.  With  such  an  approximation  made,  the  expres- 

2 

si  on  (7)  for  q reduces  to 

V 


2 2 2 2 

9!  = q_-|  - p /cA 


(16) 


2 2 

where  c.  = B /4irn  m is  the  Alfven  speed.  The  function  g defined  by  (9) 
Moo 

reduces  to 


1 

, 2 , 2.2 

ko  + p /c# 


kv  0~  + k JB 
x x y y 


(17) 


which  expression  can  then  be  used  to  calculate  the  components  of  wave 
magnetic  field. 

5.0  TRANSVERSE  COMPONENTS  OF  WAVE  MAGNETIC  FIELD 

The  combination  of  (10),  (11),  (16)  and  (17)  leads  to  the  following 
expression  for  the  transverse  magnetic  field  components  in  the  hydro  - 
magnetic  approximation: 


B 

V 


4-rcv/ C 


T 
k + 


~2 

P /c« 


k k 

O V 

T 


■*7 


+ JL 


k + p /c. 


(18) 


The  second  term  in  (18)  can  be  called  the  "bare  current"  term.  It 

2 2 2 i 

has  the  effective  plasma  dielectric  loading  of  the  factor  (k  + p /cA) ”, 
which  is  independent  of  the  angle  between  the  propagation  vector  k and  the 
DC  magnetic  field  B0.  This  term  corresponds  to  a part  of  the  magnetosonic 

'V • 


(or  fast)  wave  which,  in  cold  plasmas,  propagates  at  all  angles  to  B0  with 
the  same  phase  velocity.  For  this  reason  we  may  anticipate  that  the  contri- 
bution to  from  this  term  will  decay  as  1/r,  where  r is  the  distance  from 
the  beam.  This  contribution  therefore  is  expected  to  be  small  in  the  radia- 
tion zone. 

The  first  term  in  (18)  arises  from  the  divergence  of  the  perpendicular 
beam  current  (k  • JB).  It  has  the  dielectric  loading  factor  [fl<  +p  /c.)* 

2 2 2 l ^ M 

(k  +p  /cj]  . The  individual  factor  involving  k represents  the  transverse 
Alfven  wave,  which  is  guided  along  the  magnetic  field.  We  may  anticipate 

that  this  term  will  yield  a contribution  to  B that  does  not  decrease  rapidly 

J v 

with  distance,  at  least  parallel  to  B0.  We  shall  compute  the  contribution 

1 % 

of  this  term,  which  we  denote  as  B^,  in  the  next  section. 

5.1  CALCULATION  OF  B1 


If  one  takes  the  inverse  Fourier-Laplace  transform  of  B^,  equation  (18), 
using  also  the  expressions  for  beam  current  of  section  3.0  ^bove,  one  obtains: 


4TTvnDeV 


i k • r 

Devi  r ap  . .3.  k k e ^ ^ 

B J_  / K _pt  /■  d k o v 

/ T — e / TT  — 2 — 2 — T~~2 — 2 — r 

J 2lt1  J (2rr)  (k  +P  /Cfl)  (k^+p  /CA) 


'A7  v o 


j f 

*0  j 


o -ik- 

dre  ^ ^ [k  cosaz  + k Sin  az  ] G (jn  ,t  ) 

x y 


»r 


(19) 


^ B 

where  we  have  written  out  the  Fourier-Laplace  transform  of  J<  • explicitly. 
The  Laplace  inversion  integral  is  taken  along  the  standard  Bromwich  contour, 
i.e.  to  the  right  of  all  singularities  in  the  complex  p-plane. 

To  facilitate  the  integrations  indicated  in  (19)  we  introduce  the  cy- 
lindrical-coordinate representations  (with  B0  defining  the  z-direction) 

% 

kx  = k^  cos <p  , ky  = ^ sin<}>  \ 
x=  r cos  0 , y = r sin  0 f 

k*r  = k z 

^ o 


+ k^  r cos  (©-4>) 


(20)  . 
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k [k^cosAz  + k sinAz  ] = -=-  cos  (az  -$)  e* 


The  appearance  of  the  Heaviside  unit  step  functions  in  the  expression 
for  G (^,t)  limits  the  range  of  the  z integral  to  0 - z - Vnt  . The 
double  integration  dt  dz  can  be  integrated  by  parts  to  yield 


°°  ..  * » V r* 

f dt  e"pt  f dz  F (z  ,x  ,y  ) = J dt  e'pt  F ,y  ) ^ 


3 - 

where  F denotes  the  integrand  following  d r in  (19).  The  integration 
over  dx  dy  can  be  performed  since  the  Fourier  transform  of  the  Gaussian 


f -2  , ' 


(ri)h  f dxVlk*X  e"A  " (X  ‘RcosAVl»t  ) = 0 -V/4  e -ik^Rcosfit^ 


where  we  used  the  identity  aV||  = Q. 


-i  WiReV.V,,  r.  pt  r - nf'  C k dk  f kfdk. 

B (r,t)  = ^ dt  e’pt  J -£-£  J 

y~2  c 2irip  o -oo  2tt  o (2it) 


2 2 2^ 

e kj-A  lk°^Z  V|,t  ^ f "iv<*>  - ik  rcos(0-<)))-ik,Rcos(nt/,-<ti) 

J --*-e  cos  ( fit  -<#> ) e ■L 

,2  222  2 2 n 2tt 

(k  + P /CA)(k  + p /cA) 


where  wi  have  assumed  interchangeability  of  the  orders  of  integration  over 
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To  perform  the  ^-integration  in  (23)  we  note  that 


k^cos(fit  -4>)e 


-ik  Rcos(fit  -41)  -ik  Rcos(fit  -ij>) 

± = i — e x 

1 3R  e 


and  use  the  well  know  identity 


±ixcos0 


= Y 1 Jn  <x)« 


n = -°° 


The  ^-integration  in  (23)  can  then  be  carried  out  to  give  the  result 


v+1  -ive  i 

i e To  2-,  Jn+V  (k±r)  Jn  (k±R)  e 


n = -oo 


where  = fit  - 


The  Gegenbauer  summation  theorem  can  be  applied  to  (26) 
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Y J +n  (kx r)  Jn  (kj_R)ein  = vJ^  (kj_p) 


1 >1  1 v / 2 


where 


2 2 2 

P = r + R - 2rRcosii< 


p‘  = (r-Re"1’^)  ( r-Re1 ^) 


(28a) 


(28b) 


The  expression  for  B ((-,0  can  be  rewritten  as 


4.t:  - 


dkoko  _1V 

—n e 


oo 

/ 


where 


• 2 2 

dk±kj  e-k^  * /4 

2 u (kQ+p2/cA2)  (k2+p?kA2) 


C = V. | t-  - z 


r-Re*1 ^ v/2 


The  quantity  £ defined  in  (30)  will  be  restricted  later  to  positive 
values.  We  note  that 


. fr-Re^  V/21  _ !_  oJv(l<-tl,) 
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5.2  Backward  Wave  Radiation 

A particle  beam  preferentially  radiates  hydromagnetic  waves  propagating 
in  the  direction  opposite  to  that  of  the  beam.  Therefore,  we  confine  the 
subsequent  analysis  to  the  backward  direction,  so  that  in  (30)  z = -|z| 
and  c = V||t'+|z|>  0.  We  next  interchange  the  kQ-  and  the  kjj-  integrations 
in  (29),  and  write  the  integral  over  kQ  as 


dk 

2tt 


le4,1 

0 


dk  k e'lk°? 
o o 


(k£  + p2/C£)(k2+P2/C£J  ‘ \ 2 it 


C^7cj  k^ypvc  l 


(32) 


The  form  of  (32)  clearly  separates  the  Alfven  and  magnetosonic  contri- 
butions. We  can  evaluate  the  integral  by  Cauchy's  Theorem,  with  the  closure 
on  an  infinite  semicircle  in  the  lower  half  of  the  complex  k - plane.  This 


picks  up  residues  at  the  simple  poles  kQ=-ip/CA  and  kQ=-i (k|+p2/C|) 
producing 


1/2 


thus 


j = -2kj~[e"PE/Cfl  - e‘{(kl+  PZ/CA),/2] 


(33) 


This  then  leads  to  the  following  expression  for  B1  (r,  t) 


B1 (r,t)  _ 

v ‘ y[2  C 


_^nBeV  | I • iv+1e- 


”•  is,/1  i S'r"1"' 


3 

3R 


PJy(kj_p) 


r-Re 


1 Vtl/ 


-P5/Ca  _ e-5(k^p2/C2) 


1/2 


(34) 


If  the  spread  in  guiding  centers  a — ►O,  that  is  for  the  case  of  an 
infinitely  small  beam  cross-section,  the  Alfven  contribution  involving  the 
term  exp(-p?/CA)  in  (34)  would  lead  to  a divergent  integral  over  k^.  The 
magnetosonic  term  involving  exp  ^-s(k|+p2/Cjp  pleads  to  a convergent  k^- 
integratior,  even  in  the  limit  A — ►O. 
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5.2.1  Alfven  Contribution 

Let  us  denote  by  BlA  the  Alfven  contribution  to  (34).  The  integration 

V 

over  ki  can  be  carried  out  to  yield 


where  we  have  interchanged  the  order  of  p-  and  t1-  integrations.  By  closing 
the  contour  of  the  Laplace  inversion  integral  by  an  infinite  arc  in  the 
left-hand  p-plane,  Cauchy's  Theorem  yields  the  expression 


B1A 


( r,t)  = 


1-e 


-p2/a2 


r-Re 


1 VTp 


H(t-|z|/CA) 


(36) 


where 


vnDeV | V | | 

A a ..iJJL-tV 
c 

I z I /^a 

,+V|l'CA 


(37) 


and  we  have  set  the  (arbitrary)  phase  angle  e=0.  The  unit  step  function  in 
(36)  assures  that  no  wave  signal  arrives  at  the  point  |z|  prior  to  a time 
t = |z|/CA>  With  the  aid  of  the  identity 


. 1 ? 

[ l-e*^2] 

2e-iv*e-P2/A2  (38) 

. r-Relv*  J 

ivR  3 

. r-Reiv*  J 

A2 

then  (36)  may  be  written 


fit 


B!A(r>  *)- 


f l-e-p2/A2  l-e-(r-R>2/A21  _ 2A  f 

[ r-Re1  v^°  r-R  J ^ J 


d»pe 


-i  vt(i-p2/A2 


(39) 


where  <p  = fit  and  the  step  function  H(t-|z|/C.)  is  understood.  On  the  axis 

O O 

of  symmetry  of  the  beam,  one  has  r = 0,  p = R,  so  that  the  Alfven  wave 
magnetic  field  can  be  written: 
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kz. 


t] 

_A 

i vfiR^" 

(,.e-R2/62) 

' / . 

r = 0 

(40) 


A ( p-i v^o  \ 

i V "V 


If  R » A,  i.e.  the  mean  gyroradius  of  beam  ions  is  much  larger  than 
the  spread  in  guiding  centers  caused  by  electrostatic  expansive  effects 
on  the  originally  narrow  beam,  then  one  has  the  dominant  term 

SvA  k- 

r = o 

In  this  limit,  the  on-axis  disturbance  consists  of  a wave  contribution 
oscillating  with  frequency  w0  = n/(l  •+  V»  /C^),  and  a DC  contribution. 

At  a distance  off  the  axis  equal  to  the  mean  beam  gyroradius,  i.e. 
r = R,  one  has  ^j2  = 2R2 ( j.-cosip ) . Then  Qt 


'“l  A 

Bv  (frO 


A 


,1v*o 


r-R  ivfiR2  i-elv^° 


i-e 


- 7 (i-COS^o) 


A e"  nl- 

uK2  3n 


jd  xpe 


(41) 


ncosij/ 


(42) 


where  n=2R2/A2.  For  R»A,  the  last  term  in  (42)  will  dominate.  The  integrand 
in  this  term  is  periodic,  such  that  the  value  of  the  integral  accumulates  as 
fit0  advances  through  successive  integral  multiples  of  2u.  If  we  keep  only 
the  last  term  in  (42) 


~ 1 A 

B (t.  t) 


r=R 


•v.  2L  e"±  / dte*05* 

- n,Az  3n 


(43) 


Since  we  are  obviously  interested  in  beams  with  at  least  several 
helical  turns,  one  can  tentatively  evaluate  (43)  for  nt0=2Tr(N+e) , where 
N=l,2,3 integer,  and  0 < e < Then 


2ttN+2it6 


N 2TTm 


d^05*  + 


oirm-tZire- 

r 

di|>e 


ncosi|> 


[ d»encos*  * | 

2 TT 

1 


= N j d^e 
0 

2irfr 

= 2-n  N 1 0 ( n ) + | dipe 
0 


d^e 


ncosij> 
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where  we  have  used  the  integral  representation  of  the  modified  Bessel 
function  IQ(n).  Thus 


giA 


(r,t) 


r=R 


* _ -*NA 

- fiA2 


I (n)+  i 

1 SIT 


2«it 

— J c 


di^e 


ncosik 


For  n = R2/a2  » 1,  one  has  approximately 


Aftt0 

/ITflAR 


(44) 


At  the  fixed  position  |z|,  for  t > |z|/c^,  this  result  (44)  represents 
a pulse  whose  amplitude  grows  linearly  in  time  as  the  number  of  turns 
N = f2t/2n  of  the  beam  increases.  Also,  one  notes  that  the  amplitude  is 
larger  than  that  on-axis  (r=0)  by  the  factor  R/A>>i. 


For  r >>R,  p2/A2>>i,  so  that 


. iA 

Bv  Oo*) 


r>>R 


A 

ivfiR2 


(45) 


Thus,  the  amplitude  of  the  Alfven  pulse  outside  the  cy  linder  defined  by 

2 

r = R in  the  backward  direction  decays  as  1/r  . We  conclude  that  the 
largest  Alfven  response  occurs  at  r = R,  and  is  confined  to  a spatial 
width  about  r = R of  order  A. 


18 


5.2.2  Maanetosonic  Contribution 


To  estimate  the  magnetosonic  contribution  to  BvJ(j(;»t),  we  let  A -*■  0 in 
equation  (34)  and  carry  out  the  kj_-integration  to  obtain 


“ , s dpept  dt,  e~pt'  'p£/c^ 

Bv  (£’t}  = _AaR  W-  1 r.ReWK- 


D 2 2k 

t-(5  *f  r 

CA 


UV)H 


If  we  confine  our  attention  to  the  radiation  zone  |z|»r,  £»p,  then  to  the 
lowest  order  in  p/K  we  obtain 


A , ( [ dtl  , -■><*'  + ^A> 

2ca  3R  2tt1  5 


The  derivative  with  respect  to  R and  the  Laplace  inversion  can  be  done 
immediately  to  obtain 


- 1M  . . . A dt'  -W<Ji  r/r  \ 

Bv  (*>t)  = 2^  T~  «(t-t  -5/Ca) 


where 


®'M  ' 2h:  sn5’H(t-iz|/cA) 


V, it  + | z | 
= 1 + V,  |/c, 


Thus  we  see  that  the  magnetosonic  contribution  decays  as  l/|z|  along 
the  axis  and  is  small  compared  to  the  Alfven  term  for  |z|  >>  R. 


5.3  BARE  CURRENT  CONTRIBUTION 


We  may  use  the  same  techniques  to  evaluate  the  "bare  current"  contribu- 
tion, B (f,  t),  represented  by  the  final  term  in  equation  (18).  This  can  be 
reduced  to  the  expression 


«»%« ■ 


/ 2 v/z 

(c  +p  ) 


.V'* 


(49) 


where  we  assumed  A -*  0 since  no  divergence  problem  existed  with  the  kj_  - 
integration.  For  c » p,  the  lowest  order  contribution  to  (49)  is  given  by 


;-1^0 

Try 


H(t  - 


|z|/CA) 


(50) 


AIM 

which  is  simply  twice  the  contribution  from  B (j^.t)  given  by  equation  (48) 
above. 

5.4  SUMMARY  OF  RESULTS 

The  ion  beam  radiates  two  types  of  hydromagnetic  signals.  The  magneto- 
sonic  signal  spreads  out,  probably  almost  isotropically,  so  that  its  amplitude 
decays  roughly  as  the  inverse  distance  away  from  the  beam  source. 

The  Alfven  signal  is  confined  approximately  to  the  cylindrical  volume 
r £ R projected  by  the  beam  in  z > o into  the  backward  direction  z < o.  This 
signal  contains  both  left  and  right  hand  polarizations.  For  |r-R|>A  and 
r < R,  the  amplitude  of  the  Alfven  contribution  is  given  approximately  by 
equation  (41).  It  is  roughly  constant  across  this  projected  beam  cross-section, 
and  it  does  not  decay  with  axial  distance  |z|  away  from  the  beam  source. 

For  1 1 — R | <A,  the  amplitude  of  the  Alfven  disturbance  is  given  by  equation 
(44),  and  is  an  increasing  function  of  the  number  of  turns  made  by  the  gyrating 
ion  beam.  The  total  Alfven  signal  has  a sharp  peak  at  r=R.  This  peak  results 
from  the  anisotropic  dielectric  properties  of  the  background  plasma.  Magnetic 
field  lines  generated  by  the  gyrating  beam  current  cannot  spread  out  isotropically 


h 


T 
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into  the  surrounding  plasma  medium  as  they  would  into  a surrounding  vacuum. 
The  Alfven  signal  does  not  transport  such  field  lines  in  a direction  per- 
pendicular to  the  background  DC  magnetic  field  £ , at  least  in  the 
approximation  of  cold,  linearized  plasma  used  here;  to  the  contrary,  within 
our  approximations,  the  Alfven  signal  only  propagates  information  parallel 
to  £ . This  is  the  physical  reason  that  the  Alfven  contribution  to  the 
beam-generated  magnetic  flux  in  the  backward  direction  is  confined  to  a 
volume  |z|>o,  r<R  defined  by  the  dimension  of  the  beam  ion  gyroradius  and 
spread,  A. 
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6.0  POWER  ESTIMATES 


We  now  may  form  some  rough  estimates  of  the  total  power  radiated  by 

an  ion  beam.  We  shall  express  the  results  in  terms  of  the  total  parallel 

9 

beam  current  in  statamperes  (1  amp  = 3 x 10  statamp) 

In  = figeVn 

The  Poynting  flux  is  approximately 

S ■ B/4tt  (51 ) 


since  |E|  ^ c^|B|/c. 

6.1  Magnetosonic  Power 

For  simplicity  we  shall  assume  that  the  magnetosonic  contribution 
spreads  out  isotropically  in  the  backward  hemisphere  whose  area  we  approxi- 
mate  as  ^ 2it£  (t  ).  Then  we  may  write  the  amplitude  factor  A in  the  form 

| A [ = InVL/Kf2c  (52) 

The  time-averaged  power  follows  as 

PM-2^  'nY/c2  <53> 


6.2  Alfven  Power 


We  can  estimate  the  powers  in  the  Alfven  contribution  for  |r  - R|>A 
and  |r  - R| < A separately . For  r>R  the^power  is  negligible.  For  |r  - R|>A, 
we  estimate  the  effective  area  as  'v  irR  , and  use  (41)  to  find 

2 2 2 

P„(|r-R|>A)  'v  cAIllVJ-  = cAr“  (54) 

A 2-7~2—  T~ 

2R  n c 2c 

since  Rsi  = Vj_. 
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For  |r  - R|<  A,  we  assume  the  effective  area  9.  2ttRa,  and  use  (44) 
to  find 

2 2 

P n ( i r- R i < &)  -v  CAIn  <nto>  R (55) 

A ^ 9 A 

c 2ir 

Thus,  one  finds  the  power  in  the  Alfven  signal  given  by  (55)  exceeds 

2 

the  total  magnetosonic  power  by  a factor  ^2tt  (^tQ)  R/A  which  can  be  as  large 
as  1,000.  However,  (55)  is  nonoscillatory . 

6. 3 Power  Levels 


At  ionospheric  heights,  and  at  the  magnetic  equator  inside  the  plasma- 

7 

pause  near  3 Rc,  the  Alfven  speed  can  be  about  3 X 10  cm/sec.  If  we  express 
the  current  in  (MKS)  amperes  (In  = 3 X 10  1^  where  1^  is  amperes),  the 

power  in  both  magnetosonic  and  propagating  Alfven  waves  (for  |r  - R|>a)  is 

2 2 2 

P 9,  cAI ! i/c  = 0.03  1^  watts  (56) 

If  the  ion  beam  is  launched  near  the  geomagnetic  equator,  then  the 
Larmor  radius  of  a proton  beam  will  be  approximately  R 9 c^/Q  = 3 X 10  cm, 
using  Bq  9,  10’  Gauss.  The  amplitude  of  the  Alfven  contribution  to  the 
disturbance  field  would  be  9,  0.03  1^  in  gammas  (ly  = 10  G).  The  magneto- 
sonic wave  would  spread  out  and  have  a negligible  amplitude  by  the  time  it 
reaches  the  top  of  the  ionosphere.  The  propagating  Alfven  wave  (not  the  pulse) 
should  be  collimated  partially  by  the  magnetic  field  j3Q,  and  should  reach 
the  ionosphere  with  most  of  its  power.  As  the  beam  moves  down  the  magnetic 
field  lines  from  the  equatorial  region,  the  local  ion  gyrofrequency  changes 
due  to  increasing  geomagnetic  field  intensity.  If  a narrow  frequency  band- 
width is  desired,  only  the  spatial  region  of  about  9,  1R^  will  be  an  effective 
radiation  length.  Thus,  the  beam  should  be  turned  on  for  a time 
t < R^/c^  9.  20  sec.  The  total  energy  in  the  Alfven  contribution  under  this 
condition  will  be 


2 

W -■  P^t  9,  0.3  1^  joules 


If  the  beam  is  launched  in  the  ionosphere,  the  magnetosonic  wave  should 
couple  most  effectively  into  the  ionospheric  Alfyen  wave  guide.  The  Alfven 
wave  simply  will  propagate  along  the  magnetic  field  and  into  the  E-region. 
Since  the  beam  will  cross  the  wave  guide  region  in  about  one  second,  the  total 
energy  would  be  0.02  1^  (joules). 

Let  us  suppose  that  we  wish  to  excite  the  ionospheric  waveguide  in  such 
a way  that  a signal  strength  or  disturbance  field  6B  _ 10  mY  = 10"^  gauss  is 
present  in  that  volume.  The  total  r.m.s.  energy  present  within  the  auide 
would  be 


W = 4ttRc2  A R % 2,000  joules 

I DTT  t 

for  'v  200  km.  In  order  for  the  Alfven  wave  to  deliver  this  energy,  the 
beam  current  would  have  to  be  1^  ^ 8 amperes,  and  this  delivered  to  a potential 
of  perhaps  10  kV.  Thus  the  beam  energy  would  be  ^.0  8 to  several  megajoules, 
depending  upon  exact  conditions.  On  the  other  hand,  direct  excitation  by  the 
magnetosonic  wave  of  a similar  level  of  disturbance  would  requi re  a beam 
current  1^  ^ 240  amperes  and  a beam  energy  on  the  order  of  3 or  more  megajoules. 
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7.0  BEAM  SPREADING 

If  a well-collimated  ion  beam  at  any  of  the  current  strengths  we  have 
calculated  were  launched  from  a spacecraft,  either  at  higher  altitudes  at 
the  geomagnetic  equator,  or  within  the  ionspheric  Alfven  waveguide,  the 
initial  beam  number  density  would  surely  exceed  the  local  plasma  electron 
density  by  some  orders  of  magnitude. 

Two  possible  mechanisms  can  be  envisaged  that  would  cause  a subsequent 
spatial  spreading  of  such  an  initially  over-dense  ion  beam. 

The  first  of  these  is  space-charge.  If  the  ion  beam  is  not  neutralized 
immediately  by  either  background  plasma  electrons,  or  artificially  produced 
electrons  emitted  from  the  spacecraft,  space-charge  fields  of  the  ions  in 
the  beam  will  cause  the  beam  to  spread  on  a time  scale  of  a few  ion  plasma 
periods.  Since  the  ion  plasma  frequency  greatly  exceeds  the  ion  gyrofre- 
quency  either  in  the  upper  ionosphere  or  in  the  equatorial  magnetosphere 
at  distances  3-4  R^,  then  the  time  scale  of  the  ion  plasma  period  is  much 
shorter  than  that  of  the  Larmor  gyration  for  beam  ions.  This  implies  that 
the  spatial  spreading  of  such  beams  caused  by  space-charge  forces  will  occur 
early  in  the  first  helical  turn  of  the  gyrating  beam. 

Even  were  the  beam  to  be  charge-neutralized,  it  still  would  constitute 
an  ion  current  in  the  background  plasma.  Such  an  ion  current  would  certainly 
excite  a variety  of  beam-driven  plasma  instabilities.  Simple  calculations 
show  that  the  most  likely  instabilities  would  be  electrostatic  modes  with 
frequencies  comparable  to  the  ion  plasma  frequency  and  wave  vectors  nearly 
perpendicular  to  the  local  magnetic  field  Bo  and  lead  to  an  effective 
spread  in  guiding  centers  as  we  assumed  in  the  beam  model  of  equation  (12) 
in  section  3. 

Such  beam  spreading  probably  would  continue  until  the  local  density  of 
the  spread  beam  decreases  to  a value  roughly  comparable  to  the  density  of 
background  plasma.  A trivial  calculation  of  the  spread  a used  in  equation 
(12)  et  seq.  can  be  made  on  the  basis  of  such  a model.  For  spread  a,  the 
beam  density  (cnf^)  is 
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If  we  are  interested  in  geomagnetic  equatorial 
with  I £ '•o  80  amperes,  then  I||  a,  2.4  x 10^  e.s.u. 
cm/sec  and  Nq  102  cm”3  then  one  calculates 

4 

A ^ 5 x 10  cm 
R ^ 106  cm 
A/R  * 0.05  = 1/20 

Since  this  value  of  a/R  <<  1,  beam  spreading  appears  unlikely  to  greatly  dis- 
tort the  spatial  structure  of  the  beam  or  disrupt  significantly  the  coherent 
gyro-motion  of  the  beam.  That  is,  it  will  not  distort  the  coherence  of  Lamor 

g 

phases  of  beam  ions  necessary  to  produce  a strong  which  is  the  component 
required  to  radiate  the  waves. 

8.0  BEAM  CHARGE  AND  CURRENT  NEUTRALIZATION 

8.1  Charge  Neutralization 

It  is  clear  a prior  that  for  beam  currents  required  in  the  previous  sec- 
tions, that  is  currents  of  ions  on  the  order  of  tens  of  amperes,  there  will 
be  severe  charge-neutralization  problems  encountered  by  a spacecraft  plat- 
form used  to  launch  such  beams. 

The  beam  voltages  required  to  produce  beam  speeds  V j | of  several  times 
the  local  Alfven  speed  cft  are  tens  of  kilovolts.  Unless  negative  charge 
equal  to  that  of  the  ion  beam,  i.e. 

t 

Q.-f  Iudt 

0 

is  emitted,  a spacecraft  potential  of  many  kilovolts  would  ensue. 

The  local  plasma  at  the  geomagnetic  equator  near  R ^ 3 to  4 R^.  cannot 
supply  sufficient  ion  return  current  to  neutralize  the  spacecraft  charge. 

For  this  reason,  one  would  require  electron  emitters  on  the  spacecraft  to 
supply  a sufficient  negative  current  to  equalize  the  negative  charge  caused 
by  positive  ions  leaving  the  spacecraft. 


launching  of  the  beam 
If  we  have  V( | ^ 2C^  ^ 103 


This  magnitude  of  electron  current  (I  % 80  amps)  can  be  supplied  at 

relatively  low  voltages  from  large  cathodes.  For  example  large  Barium- 

3 2 

Strontium-Oxide  coated  cathodes  of  dimensions  ^ 3 x 10  cm  have  been  con- 
structed at  TRW  for  use  in  laboratory  plasma  chambers.  These  cathodes  are 
easily  capable  of  producing  electron  emission  currents  of  hundreds  of  amperes 
at  50-100  volts  on  the  extractor  grid. 

Thus,  charge-neutralization  of  the  spacecraft  appears  to  be  within  state- 
of-the-art  technology,  but  it  significantly  adds  to  the  experimental  energy 
budget.  80  amperes  of  electron  current  at  -v  100  volts  consumes,  even  with 
one-hundred  percent  efficiency,  about  8 kW.  This  is  comparable  to  the  power 
in  the  ion  beam  itself.  In  addition,  one  would  have  to  include  the  heater 
power  for  this  electron  gun,  which  will  be  several  kW. 

8.2  Current  Neutralization 

It  was  suggested  in  section  3.0  above  that  beam  current  neutrality 
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(Jjl  = 0)  can  be  achieved  using  the  same  source  of  electrons  as  that  used 
to  charge-neutralize  the  spacecraft  itself. 

The  picture  painted  there  is  one  in  which  the  ion  beam  is  injected  at, 
say,  an  initial  pitch  angle  of  about  45  degrees  to  J3o  , and  thus  describes 
a helical  current  element  of  mean  radius  R ^ 10^-  10 & cm.  The  electron 
emitter  furnishing  the  neutralizing  current  is  then  assumed  to  inject  its 
electrons  along  magnetic  field  lines  within  this  cylindrical  shell. 

However,  exact  spatial  current  neutralization  will  require  not  only 
that  ion  and  electron  beam  currents  be  equal,  but  also  that  beam  velocities 
parallel  to  B be  equal.  For  10  keV  proton  beams,  one  then  would  require 

<v,° 

'v  5 eV  electron  beams  to  match  velocity.  With  reasonably  sized  cathodes, 
drawing  80  amperes  of  electrons  at  only  5 eV  is  limited  by  beam  perveance 
unless  large  beam  cross-sections  are  employed.  Such  technological  pro- 
blems are  not  considered  to  be  beyond  the  realm  of  solutions. 

Insofar  as  neutralization  by  ambient  plasma  reaction  currents  is  con- 
cerned, our  formulation  of  the  original  problem  via  equations  (1)  et  seq . 

p 

contains  the  quantity  J , the  plasma  reaction  current  density.  This  self- 
consistent  quantity  describes  the  plasma  reaction  through  the  plasma  dielec- 
tric function,  and  gives  the  linearized  shielding  of  the  beam  charge  and 
current  due  to  the  ambient  plasma. 
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9.0  BEAM  MODULATION 


Energy  requirements  for  the  ion-beam  ULF  generator  operated  from  a 
spacecraft  located  around  the  geomagnetic  equator  at  some  3-4  R£  are  est- 
imated at  several  megajoules,  for  operation  of  neutralized  beams  at 
some  80  amperes  and  10  kilovolts. 

It  the  beam  is  launched  in  a region  where  B0  ^ 200-y,  then  £}/2tt  ^ 3 Hz 
and  one  turn  of  the  beam  takes  ^ 0.3  sec.  For  ten  tums>  one  requires  3-second 
operation,  so  that  beam  energy  would  be  'v  2.4  megajoules. 

Modulation  of  such  a beam  would  not  lead  to  any  advantage  because  of 
its  necessarily  short  duration. 

10.0  HEAVIER  ION  BEAMS 


Again,  because  the  beam  energy  required  is  fixed  by  the  Alfven  speed 
in  the  ambient  (assumed  proton)  plasma,  one  cannot  justify  using  heavier 
ions  to  form  the  beam.  We  have  found  that  VM/cA  >2  is  required,  and  thus 
to  produce  this  condition  using  beam  ions  of  atomic  number  A simply  mandates 
an  increase  in  beam  energy  A-times  greater  than  that  for  a proton  beam. 

For  this  reason,  use  of  heavier  ion  beams  is  considered  to  be  a distinct 
disadvantage. 

One  might  argue,  on  the  other  hand,  that  the  use  of  a heavier  ion  such 
as  Cesium,  with  a much  lower  ionization  potential,  could  lead  to  an  offsetting 
gain  in  energy  saved  in  producing  the  ions.  However,  given  even  equal  ionizat- 
ion efficiencies,  the  lightest  ion  of  low  ionization  potential  is  Lithium  at 
^ 5.4  eV,  compared  with  Hydrogen  at  13.6  eV.  This  ratio  is  13.6/5.4  = 2.5 
while  A=7  for  Li.  As  A increases  the  disparity  between  the  ratio  of  ionization 
potentials  to  the  ratio  of  required  beam  energies  to  produce  the  specific  vel- 
ocity Vi i >2  c A increases. 


Energetics  alone  militate  against  the  use  of  heavier  ion  beams. 
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11.0  CONCLUSIONS 


The  following  general  conclusions  are  based  on  the  results  of  this 
investigation  of  the  efficacy  of  generating  ULF  radiation  from  ion  beams 
emitted  from  satellites,  and  based  on  those  conclusions,  the  recommendations 
concerning  further  pursuit  of  this  concept  are  enumerated  in  Section  12. 

1.  The  calculation  of  hydromagneti c radiation  from  a proton  beam 
injected  into  a dilute,  magnetized  background  plasma  has  been 
carried  out  in  linearized  analysis. 

2.  For  such  beams  emitted  near  the  geomagnetic  equator,  the  space- 
craft must  reside  within  the  plasmasphere  so  that  sufficient 
ambient  plasma  is  present  to  allow  coupling  of  beam-generated 
ULF  into  the  backward  plasma  wave  at  frequencies  u>  < n = 

beam  gyro-frequency. 

3.  Backward-propagating  disturbances  in  the  Alfven  mode  have  been 
found.  These  have  the  property  of  being  largely  confined  to  a 
cylindrical  volume  r <_  R defined  by  the  mean  gyroradius  of  the 
ion  beam  (emitted  into  z<  0). 

4.  Almost  isotropically  propagating  magnetosonic  wave  disturbances 
are  also  found. 

5.  For  a given  beam  current  (or  power)  the  power  in  the  guided 
backward  Alfven  disturbance  can  exceed  the  power  in  the  iso- 
tropic magnetosonic  disturbance  by  a factor  on  the  order  of 

1 ,000. 

6.  Given  a unity  efficiency  of  guidance  of  the  beam-generated 
Alfven  disturbance  along  geomagnetic  field  lines  down  to  the 
ionosphere,  we  calculate  that  to  produce  a ULF  disturbance  6B 
of  10  my  throughout  the  earth-ionosphere  waveguide  requires  an 
r.m.s.  energy  of  2 kilo-joules  in  the  guide,  and  this  translates 
into  a beam  energy  of  the  order  of  using  an  efficiency  factor 

2 x 10'3. 


7.  Typical  operating  parameters  for  such  beams  are  Ig'v-  29  amperes, 

VgVIO  k V,  Pg  * 0.2  MW,  Wg  * 2 MJ. 

8.  A neutralizing  electron  emitter  is  required  to  prevent  the  space- 
craft from  attaining  large  negative  potential.  If  a low-voltage, 
high-current  electron  emitter  were  used,  say  I ^ 20  amperes  at 
Vg  = 5V,  only  100  watts  of  electron  beam  power  would  be  expended. 

Such  emitters  appear  to  be  within  the  state-of-the-art.  Heater 
power  consumption  would  have  to  be  added  to  the  power  budget,  at 
an  expense  of  about  10  eV  per  electron. 

9.  Electrostatic  effects  can  be  expected  to  cause  an  initially  narrow, 
intense  collimated  ion  beam  of  small  cross-section  to  expand  radially 
about  its  mean  radius  of  gyration,  R = V /fi.  The  spread  A so  produced 
is  estimated  on  the  assumption  that  the  beam  will  expand  until  its 
number  density  is  roughly  equal  to  the  background  plasma  density,  and 
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that  for  the  beams  considered  here,  A/R  ^ 0.05,  where  R^3-10X10  cm. 
This  amount  of  expansion  would  have  negligible  effects  on  the  rad- 
iation due  to  the  perpendicular  component  of  beam  current  density, 
J®. 

10.  Plasma  reaction  currents  are  included  in  our  calculation.  They  cont- 
tribute  the  plasma  dielectric  loading  that  controls  the  hydromagneti c 
radiative  characteristics  of  the  beam. 


11.  Because  of  the  energy  storage  requirements  calculated,  it  is  infeasible 
to  consider  phase,  velocity,  or  density  modulation  of  such  an  ion  beam. 

12.  Also  because  of  energy  considerations,  it  appears  to  be  distinctly 
disadvantageous  to  use  heavier  ion  beams.  The  beam  velocity  is  fixed 
by  the  magnitude  of  the  A1 fven  speed  in  the  ambient  plasma.  Thus,  to 
obtain  the  necessary  beam  velocity,  the  lowest  energy  requirement  is 
for  a proton  beam. 

13.  In  any  practical  consideration  of  total  energy  storage  requirements 
for  such  a spacecraft-borne  ion  beam,  one  must  include  the  ionization 
energy,  dissipation  in  the  power  supply  used  to  apply  accelerating 
voltage,  and  energy  to  produce  the  neutralizing  electron  current.  The 
ion  beam  energy  alone  is  estimated  to  be  ^ 2MJ.  Thus,  one  could 
anticipate  a 3MJ  storage  requirement.  At  even  0.5  lb  per  kilojoule, 
capacitive  storage  would  be  on  the  order  of  1,500  lb.  Furthermore,  even 
at  a recharging  rate  of  lOkW,  which  represents  a solar  panel  of  pro- 
hibitive size,  the  recharging  time  would  be  on  the  order  of  5 minutes. 
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14.  The  final  conclusion,  based  on  this  information,  is  that  the 
production  of  ULF  signals  in  the  frequency  range  below  a few 
Hertz  by  means  of  an  ion  beam  emitted  by  a spacecraft  near 
the  geomagnetic  equator  at  altitudes  ^ 3 to  4 R£  has  very 
severe  requirements,  and  at  best  would  supDly  some  sort  of 
"bell-ringer"  signal.  However,  the  signature  of  the  signals 
we  have  calculated  would  be  difficult  to  control  under  varying 
magnetospheric  and  ionospheric  conditions. 

12.0  RECOMMENDATIONS 

Based  on  the  foregoing  calculations  and  conclusions  which  we  have 
drawn,  the  following  recommendations  can  be  made: 

1.  The  concept  originally  envisaged  of  an  ion  beam  emitted  from  a 
satellite  near  the  geomagnetic  equator  at  distances  of  some  3 
to  4 earth-radii  should  probably  not  be  pursued  further.  The 
results,  especially  the  energy  storage  requirements,  appear  to 
argue  against  any  possible  large-scale  implementation  of  such 
a system.  To  provide  coverage  would  require  several  spaced 
satellites. 

2.  Even  under  ideal  quiescent  magnetospheric  or  ionospheric  con- 
ditions, the  signals  produced  would  be  marginally  useful.  Given 
a dynamic  magnetosphere-ionosphere  system,  this  concept  appears 
most  difficult  to  utilize  in  a controlled  fashion. 

3.  Some  follow-on  studies  could  be  defined  to  examine  related 
concepts,  and  to  investigate  certain  features  of  ion  beam 
dynamics  that  cannot  be  calculated  with  any  degree  of  confidence. 
Among  these  are: 

• A study  of  possible  location  of  the  ion  beam  within  the 
ionosphere  waveguide.  Since  much  of  the  power  produced  by 
such  a beam  appears  in  the  magnetosonic  mode  directly,  and 
since  this  is  the  mode  that  one  assumes  produces  the 
leakage  field  into  the  earth-ionosphere  cavity  that  appears 
at  the  surface,  then  a calculation  placing  the  beam  source 
in  the  waveguide  is  perhaps  indicated.  This  could  best  be 
done  using  numerical  analysis  and  a model  of  the  waveguide. 
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• A laboratory  experiment  can  be  performed  in  a large  plasma 
facility  at  TRW  to  study  some  features  of  ion  beam  dynamics. 
However,  one  cannot  scale  the  ULF  radiation  because  the 
hydromagnetic  wavelengths  exceed  the  dimensions  of  laboratory 
chambers.  Hence,  small  scale  beam  dynamics,  such  as  electro- 
static stability,  expansion,  etc.,  would  be  the  objectives 

of  such  a laboratory  investigation. 

t Numerical  modeling  of  the  beam-plasma  system,  using  codes 
written  for  that  purpose,  would  be  an  expensive  but  perhaps 
more  illuminating  way  to  further  the  analysis. 

• Creation  of  hydromagnetic  disturbances  by  MPD  arcs  should 
also  be  examined.  However,  one  can  size  this  or  any  similar 
concept  by  some  energy  considerations.  One  must  always  come 
out  with  megajoule  power  supply  requirements  in  all  such 
systems. 


